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Abstract 

We present real-space renormalization group (RG) calculations of the crit- 
ical properties of the random-field Ising model on a cubic lattice in three 
dimensions. We calculate the RG flows in a two-parameter truncation of the 
Hamiltonian space. As predicted, the transition at finite randomness is con- 
trolled by a zero temperature, disordered critical fixed point, and we exhibit 
the universal crossover trajectory from the pure Ising critical point. We ex- 
tract scaling fields and critical exponents, and study the distribution of barrier 
heights between states clS cl function of length scale. 
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The critical properties of the random-field Ising model have been the subject of 
intense controversy. One of the simplest of disordered systems, the model is governed by 
the Hamiltonian 



n = -jY.SiS,-Y.{H + K)s,, (1) 

where Sj = ±1 is an Ising spin on a cubic lattice, if is a homogeneous external field and 
the variables hi are independent Gaussianly-distributed random fields of mean zero and 
variance cr^. The notation {ij) indicates a sum over nearest-neighbor sites. The contro- 
versy surrounded the existence of a ferromagnetic phase transition at weak disorder in three 
dimensions. Supersymmetry techniques and experiments by Hagen and co-workers 
appeared to indicate that arbitrarily weak disorder would break the system up into random 
domains at long enough length scales. Simpler arguments given by Grinstein and Ma 0, 
along with experiments at Santa Barbara suggested that ferromagnetism persisted until 
a critical value of the disorder was reached. The issue was settled definitively by Imbrie and 
others 0,^, who proved that the latter viewpoint was the correct one. 

It was realized that the unusual experimental problems posed by the random-field model 
are the result of extremely slow, glassy dynamics in the system. Initial treatments of the 
dynamics P] concentrated on the non-equilibrium coarsening of the domains in the ferro- 
magnetic phase: the domain walls pin on the random fields, and as the domains grow these 
pinning barriers grow too. When the barriers become large compared to the temperature, 
thermal activation becomes ineffective and the system falls out of equilibrium without es- 
tablishing long-range order. Later analyses by Bray and Moore |10| and by Fisher |Tl 



indicated that this glassy behavior extends also to the equilibrium dynamics near the crit- 
ical point. They argue that as the temperature approaches the ferromagnetic transition 
temperature t = (T — Tc(cr))/Tc(cr) and the correlation length ^ ~ t^'^ grows, the 
effective coarse-grained Hamiltonian fiows to a zero temperature, disordered critical fixed 
point. Unlike most critical points, the equilibrium energy scale E ~ t~^'^ at the correlation 
length becomes much larger than the temperature (violation of hyperscaling) . As a result. 



the temperature is an irrelevant variable in the renormalization group (RG) treatment of 
the system. The exponent ip determining the divergence of the energy barriers B ~ t~^'^ is 
assumed to equal 6, which follows from the irrelevant eigenvalue governing the RG flow of 
the temperature. Instead of a competition between bond energy and thermal fluctuations, 
the battle at long length scales is between the bonds and the (renormalized) disorder. The 
dynamics, which proceeds by thermal activation with time constant r ~ exp{B/kT), slows 
down exponentially already above Tc{a). 

We present here a real-space RG calculation which directly confirms this theoretical pic- 
ture. Figure 1 previews our results: the pure Ising critical point Tc is unstable to random- 
field disorder, and the rescaled effective temperature and disorder flow to a disordered critical 
fixed point at T = 0, cr = ac- As well as providing a useful qualitative picture of the RG 
flows in the model, our method gives values for the critical exponents at the two fixed points. 
In addition to the simple two-parameter truncation of Hamiltonian space, we consider ad- 
ditional forms of disorder to (a) confirm that they correspond to irrelevant operators, and 
(b) improve our critical exponents. The calculation gives more precise exponents than direct 
methods, though our systematic errors are potentially large because of our small system size. 

We investigate the model by a straightforward computational implementation of the 



real-space RG [[L4 -[14| on a cubic lattice. For given values of the bond strength J, the 



external field H, and the variance of the random fields, the procedure is as follows: 

1. Choose random values for the variables hi on each site of the lattice. 

2. Divide the lattice into cubic blocks of eight spins each. For each block, we define a 
coarse-grained spin variable. 

3. For each configuration of the coarse-grained spins, sum all the Boltzmann factors for 
configurations of the original lattice consistent with those spins, using the so-called 
'majority rule' (see Ref. fl^). 

4. This defines a renormalized Hamiltonian Ti' which can be inverted to give a new set 
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of bond strengths J^j and fields h'^ on the blocked lattice, as well as a number of 
longer-range interactions between two or more spins. 

5. In the simplest case we discard all the longer-range interactions and define the renor- 
malized bond strength to be the mean J' = (Jlj) of the new bonds, the renormalized 
external field to be the mean H' = {h'^ of the new fields, and the renormalized variance 
to be the variance of the new fields cr'^ = {h'-^) — {h[)'^. Later, we will consider more 
sophisticated versions of the RG which include higher moments in the distribution of 
bonds and fields than just the mean and variance considered here. These all turn out 
to be irrelevant operators, but their inclusion in the calculation can improve the results 
for the critical exponents. 

To achieve good statistics we average the values of J' and a'^ over many different realizations 
of the randomness. For consistency with earlier work, we quote our results in terms of the 
ratios T/J and H/J, and the standard deviation of the random fields or 'randomness' a/ J. 

Employing the method first in two dimensions for a 4 x 4 system (with rules analogous to 
the three-dimensional ones described above) we find no non-trivial fixed points other than 



the pure Ising critical point |T5|. We conclude that, within this approximation, there is 
no phase transition at finite randomness in two-dimensions. This is in agreement with the 
findings of Imbrie 0, Berker |16|, Bricmont and Kupiainen and others. 



In three dimensions, the size of the system we can study is limited by step 3 above, in 
which we are required to sum over all spin configurations of the lattice. We are working 
on more sophisticated algorithms to speed the calculation, but for the moment we present 
results for a system of 2 x 2 x 4 spins. This is large enough to give reasonable results, but 
small enough for the method to be directly applicable. Our numerical results for the RG 



flows. Figure 1, are consistent with the flow diagram postulated by Bray and Moore []T0 
Each point in Figure 1 is an average over 100 different realizations of the random fields. 
In the regions close to the fixed points, we also performed a number of runs in which we 



averaged over as many as 10^ different realizations of the randomness in order to improve 



the accuracy of our values for the critical exponents. 

There are two non-trivial fixed points. One is the normal Ising fixed point at finite tem- 
perature and zero randomness, and the other is at zero temperature and finite randomness. 
This latter fixed point governs the phase transition between the paramagnetic (disordered) 
state and the ferromagnetic one. As predicted by Fishman and Aharony |T^, all three 
parameters in the problem — temperature, external field, and randomness — are relevant at 
the pure Ising critical fixed point. At the disordered critical fixed point only the external 
field and the randomness are relevant; the temperature is irrelevant. There exists a unique 
trajectory which leads from one critical fixed point to the other. This curve determines the 
crossover behavior for weakly disordered systems: far from Tc{a) for small a the system will 
have critical fluctuations given by the pure fixed point, and the growing influence of disorder 
is of a universal form given by this trajectory. 

The RG transformation does not actually provide us with a single renormalized value for 
the bond strength J', but with a distribution of strengths, which is considerably skewed from 
a pure Gaussian. We can map the RG flows more accurately if we allow for this distribution, 
parameterizing it by three quantities: the mean bond strength J, the standard deviation aj, 
and the skewness 7j. Again, the trajectories leading away from the fixed point are universal. 

While the entire trajectory contains universal information about the critical behavior, 
the most commonly measured properties are the positions of the critical fixed points, and 
the critical exponents. Our value for the critical temperature at the pure Ising critical point 
is Tc = 5.38, which is substantially higher than the known value of 4.51, due to the small 
system size. Near this point, the three parameters (T — Tc)/J, H/ J, and a/ J are expected 
to scale independently under the RG transformation, and our calculation confirms this. Our 
value for the critical randomness at the disordered critical fixed point is cXc = 1.675 ± 0.002. 
It has been conjectured that near this fixed point the correct eigenvectors of the RG 
transformation are not the bare parameters of the problem, but are instead T/J, H/J, and 
the linear combination (cr — a^/ J + A{T/J), where A is a mixing constant. We believe 
however that the T = plane should be invariant under the RG transformation, and this 

5 



would imply that A = in general. Our calculation confirms this. On the other hand, we 
see no physical reason why T/J should be an eigenvector at the disordered critical point, 
and in that case the most general scaling fields for the problem would be H/ J, {a — o"c) / J, 
and T/J + B{a — (Jc)/ J- Within our numerical calculation we find that 5 = 0, but not 
because of any symmetry of the Hamiltonian. B is zero because the fiow is asymptotically 
dominated by the lowest-lying states consistent with each configuration of coarse-grained 
spins. The argument that leads us to this result becomes invalid in the thermodynamic 
limit, so it is not clear from our calculation whether B will be zero for an infinite system. 

In the more sophisticated calculation, where we choose the bond strengths at random 
from a skewed Gaussian distribution, we find the disordered critical fixed point at T = 0, 
a = 1.55, aj = 0.36, 7j = -0.20. 

Table | shows our values for the eigenvalue exponents given by linearizing the RG fiows 
around the fixed points, together with those of previous simulations. Getting from the 
eigenvalue exponents to real critical exponents is straightforward: for example, the barriers 
are thought to diverge with exponent ipu = yj/i/a- We have included in the table estimates 
of the statistical errors. In addition to these statistical errors, we expect large systematic 
errors because of finite-size effects on the small lattice. Even with these systematic errors, 
however, our results at the important disordered fixed point are competitive with previous 
results. We are extending our work to a 4 x 4 x 4 system, which should allow us to include 
longer-range renormalized interactions and extract reliable, accurate exponents. 

The eigenvalue exponents are expected to satisfy a number of different inequali- 



ties [|n],|T8|-||]: Vh -yj < d/2, vh < d, yj < d/2, yj < d - 1, and yn - yj/2 > d/2. 
Our calculated values satisfy all of these inequalities except for the last one, where yn is 
slightly less than the required value of ^{d + yj) = 2.00. This inequality is an expression 
of the fact that the magnetization fluctuations on the phase boundary between para- and 
ferromagnetic phases are expected to diverge in the limit q —>■ [^. On the present small 
lattice it is not surprising then that this inequality fails. We hope that a calculation on a 
larger lattice would bring the value of yn up. 
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For the renormalized system with only two coarse-grained spins, it is a straightforward 
matter to calculate the distribution of barrier heights between the four possible states of 
the system. Though the small size of the lattice probably means that the exact shape 
of the distribution is not representative of the distribution in bulk systems, the scaling 
behavior of the distribution with increasing length scale, which we can extract from our RG, 
should be qualitatively correct. We calculate the barrier height distribution for a system 
with initial values J, a of the bond strength and randomness, and also for systems along 
the RG trajectory that starts at this point. The successive distributions correspond to the 
distribution for the original system on length scales increasing by a factor of, in this case, two 
for each step along the trajectory. In the inset to Figure 2 we have plotted the barrier height 
distributions for successive steps along one RG trajectory. This trajectory, which is indicated 
by the thicker line in Figure 1, starts fractionally below near the cr = axis, moves along 
the edge of the phase boundary and lingers briefly near the disordered critical flxed point, 
before turning away towards the ferromagnetic flxed point at T = 0, cr = 0. The barriers 
are plotted in multiples of the bond strength at each stage, and the vicinities of the flxed 
points are visible in the progression to longer length scales as stationary regions in which the 
distribution changes little. Figure 2 shows the mean barrier height, this time with the factors 
of J included, as a function of length scale. Scaling near the flxed points is represented by 
the straight-hne portions of this plot. The zero slope follows from hyperscahng; the slope 2 
follows from the area of an interface. The structure at shorter length scales will probably 
not have an effect on experiments, since it represents mean barrier heights much less than 
kT . But the behavior at longer scales, particularly the crossover between the disordered and 
ferromagnetic flxed points might be measurable in an experiment that correlates the size of 
magnetization fluctuations (related to the length scale) with the time-scale on which they 
occur (related to the barrier height). 

In summary, we provide a direct implementation of the RG for the random-fleld Ising 
model on a cubic lattice, a model whose unusual scaling behavior was the cause of substan- 
tial controversy and which now provides the best understood example of glassy dynamics. 
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Our results for the RG flows under coarse graining and for the position and nature of the 
fixed points confirm earlier conjectures. It is encouraging that the statistical errors on our 
exponents are much smaller than the errors generated by other techniques. If the systematic 
errors that arise from working on the smallest possible lattice can be reduced by going to 
larger lattices, the technique promises exponents of greater accuracy than those available at 
present. The RG also gives us a simple method for calculating the barrier height distribution 
as a function of length scale. The results may have experimentally measurable consequences. 

We would hke to thank J. F. Marko for helpful suggestions. This work was partly funded 
by the Hertz Foundation (BWR), the Science and Engineering Research Council of Great 
Britain (MEJN), and the NSF under grants DMR-91-18065 (BWR, JPS) and DMR-91- 
21654 (GTB). 
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TABLES 



Fixed point 


exponent 


value 


best result 


Finite T 


yr 


1.206 


1.594 ± 0.004 


fixed point 


VH 


2.212 


2.488 ± 0.004 




Va 


0.509 ± 0.002 




r = 


Va 


0.672 ± 0.005 


0.9 ±0.15 


fixed point 


VH 


1.88 ±0.01 


2.95 ± 0.05 




yj 


1.00 ± 0.05 


1.5 ±0.2 



TABLE L Results for the six eigenvalue exponents at the two critical fixed points. The figures 



in the last column are taken from Ferrenberg and Landau and Nattermann and Villain 0. 
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FIGURES 

FIG. 1. Lines of RG flow through Hamiltonian space for the case of zero external field. The 

range (0, oo) has been mapped onto the unit square by taking liypcrbolic tangents. The pure Ising 
fixed point is unstable to randomness: the universal critical behavior along the disordered phase 
boundary is determined by the fixed point at (T = 0,(7 = ctc). The trajectory marked with a 
thicker line is the one used in the calculation of the barrier height distribution, Figure 2. 

FIG. 2. The mean barrier height as a function of length scale for a system with small random- 
ness, slightly below the critical temperature. The dotted lines indicate the expected scaling of the 
mean barrier height in the vicinity of the three fixed points: the slopes are zero, 1.00 it 0.05, and 
two, respectively. The inset shows the actual distributions of the barrier heights in multiples of J 
for each point along the RG trajectory. 
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